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The specific heats of exactly solvable alternating layered planar Ising models with 

strips of width nii lattice spacings and "strong" couplings Ji sandwiched between 

strips of width 7712 and "weak" coupling J2 , have been studied numerically to inves- 

Cd tigate the effects of connectivity and proximity. We find that the enhancements of 

I the specific heats of the strong layers and of the overall or 'bulk' critical tempera- 

Ch ture, TJJi, J2]it^i,^tT'2), arising from the collective effects refiect the observations of 

O 

I ^1 Gasparini and coworkers in experiments on confined superfiuid helium. Explicitly, 

T-H we demonstrate that finite-size scaling holds in the vicinity of the upper limiting 

> 
00 critical point Tic (of Ji/ks) and close to the corresponding lower critical limit T2c 



(oc J2/kB) when mi and 7712 increase. However, the residual enhancement, defined 



t:]" via appropriate subtractions of leading contributions from the total specific heat, is 

o 

C^ dominated (away from Tic) by a decay factor l/(r7ii + 70,2) arising from the seams 

j>! (or boundaries) separating the strips; close to Tic the decay is slower by a factor 

^ ln(777l/777o). 
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I. INTRODUCTION 

Many experiments performed on ^He at the superfuid transition in various spa- 
tial dimensions^, reveal excellent agreement with general finite-size scaling theory^Sl. 
Furthermore, when small boxes or "quantum dots" of helium were coupled through 
a thin helium film, the effects of connectivity and proximity were discovered and 
quantified^i^. 

To gain some more detailed theoretical insights into the proximity effects, we 
study here the specific heats of an alternating layered planar Ising model, which 
consists of infinite strips of width mi lattice spacings in which the coupling energy 
between the nearest-neighbor Ising spins is Ji, separated by other infinite strips of 
width m2 whose coupling J2 is weaker. This is illustrated in Fig.l. 
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FIG. 1. The planar, square lattice alternating layered Ising model considered. The widths 
mi and 771-2 are measured in nearest neighbor lattice spacings while nearest-neighbor Ising 
spins (Ti = ±1 are coupled via pair Hamiltonians Jijaiaj. 



When J2 vanishes, the model becomes a system of noninteracting infinite strips 
of finite width, each of which essentially behaves as a one dimensional Ising model. 
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FIG. 2. (a) The specific heats per site of layered systems with relative strength r = J2I J\ = 
0.3 and relative separation s = 1 for m\ = 7x12 = 2, 4, 6, 8, 12 and 16 . The amplitude A{r^ s) 
of the logarithmic divergence at T^ir^ s), which dominates for m\ = 2, decreases rapidly as 
7TT-1 increases. Thus, the small spike at the "true" bulk critical point, Tc, becomes barely 
visible for rrii > 8. However, as demonstrated in Fig. 2(b), it remains evident at larger 
values of mi when r increases towards unity. On the other hand as A{r, s) becomes small 
two quite distinct rounded peaks appear and move toward the limiting values Tic and T2c 
as mi = 7712 -^ 00. 
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FIG. 2. (b) Plots of the specific heats for fixed mi = m2 = 16 and so s = 1 but for 
increasing relative strength r. Note that the logarithmic peak at the overall or bulk critical 
point, Tc{r, s), indicated by short vertical dotted lines, remains clearly visible when r = 0.7 
and still dominates entirely when r = 0.9. 



This means, in particular, that the specific heat is not divergent but rather has a fully 
analytic rounded peak. However, as long as J2 7^ 0, the system is a two-dimensional 
bulk Ising model, whose specific heat per site diverges logarithmically at a unique 
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bulk critical temperature Tc{Ji, J2', frii, 1112) in the form 

C{T)/kB oc -A(r, s) In |1 - T/T,| + ■ ■ ■ , (1) 

where we have introduced the basic weakness or coupling ratio, r, and the relative 
separation distance, s, namely 

r = J2/J1 < 1, s = m2/mi. (2) 

In fact, as will be showrP, the amplitude A{r^ s) of the logarithmic divergence 
decays exponentially as a function of the minimum of mi and 1712; indeed, at fixed 
s and r ^> 1, the amplitude decays as Qmic^*^™"^, where Q oc (1 — r)s/{l + s) as 
mi —7- 00. This behavior is evident for r = 0.3 in Fig. 2(a), which shows that the 
divergence, while obvious and dominant for nii = 1112 < 3, rapidly becomes no more 
than a minuscule spike, which soon becomes invisible on any graphical plot. On 
the other hand, for greater values of the coupling ratio r the logarithmic divergence 
remains dominant for larger values of mi and 1712 as seen in Fig. 2(b). But returning 
to Fig. 2(a) with r = 0.3, one observes that as soon as the strip widths, mi = m2, 
exceed three lattice spacings, there appear two further specific heat peaks, albeit 
rounded; these grow rapidly in height and sharpness and as mi and m2 grow in 
magnitude, they soon dominate the plots. 

Now Fig. 2 is based on exact analytic calculations expounded in Part 11 of this 
article^. In fact, the analysis of the finite-size behavior of planar Ising models based 
on the exact solution of Onsager, as extended by Kaufman^^, goes back to the work 
of Fisher and FerdinancP'^ in 1969. Specifically, the solubility of arbitrarily layered 
planar Ising models was first noted and reported at a conference in JaparP', while, 
independently, McCoy and Wil^ developed and analyzed randomly layered Ising 
models. The thermodynamics for regularly layered models was developed by Au- 
Yang and McCoy^ and Hammf^, while the scaling behavior of a single strip of finite 
width was elucidated by Au-Yang and Fisher^. 
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In general the bulk critical temperature can be simply stated, for a layered dis- 
tribution a#^ 

kBT,{{\ncoi\i{JjkBT,))) = 2{{Jy)), (3) 

where the brackets ((■)) denote an average over the distribution, random or regular 
of the distinct number (say n < oo) of lattice spacings constituting a layer of finite 
width. For the alternating layered Ising model, this becomes 

2Ji(l + rs) = kBTc[lncoth{Ji/kBTc) + slncoth(rJi/A;BTj], (4) 

which depends only on the weakness ratio r and the relative separation s. 

Then as mi and m2 become large, the upper and lower rounded peaks approach 
limiting values. Tic and T2c (as evident in Fig. 2(a)), which, in fact, match the corre- 
sponding bulk (i.e., uniform) two-dimensional Ising models with coupling constants 
Ji and J2. Thus the limiting values Tic and T2c(r) are known!^ ^^ * ^^* and given by 

ksTic/Ji ^ 2.269185312, kBT2c/Ji ^ r ■ 2.269185312. (5) 

It proves easy to establish the expected inequalities 

T2c(r) < T,{r, s) < Tic. (6) 

II. QUALITATIVE OBSERVATIONS 

To explore further and develop the analogies with the observations on superfiuid 
helium systems, we retain the value of the weakness ratio r = 0.3 (used in Fig. 2(a)) 
but increase the relative layer separation to s = 2. The results for mi = 8 and 16 
(as used in Fig.2(a)) are presented in Fig.3: see the solid curves. As anticipated, no 
sign of any singularity at T^ is visible. It should be noted, nonetheless, that were 
one to examine the overall spontaneous magnetization, Mq{T), one would find — and 



on a plot "see" — that Mq vanished identically for T > Tc but was nonzero (and 

varying as oc (T^ — T)^ with /3 = | for 2D Ising layerd^'^'^ as soon as T < Tg. In the 

experiments on superfluids the analogous statement concerns the overall superfluid 

density ps{TY] this vanishes identically above the overall or bulk lambda transition 

at Ta(= Tc) but is detectable, via setting up persistent superflow fluid currents, below 

Tx'^. (In a bulk 3D superfluid PsiT) varies as {T\ — T^ with C, — 0.67, but in a 

planar 2D superfluid, psiT) increases discontinuously at the corresponding lambda 

point, Tx, on lowering the temperature^.) 

On the other hand, the temperatures of the upper and lower rounded maxima 

increase (and decrease, respectively), as mi increases in Fig. 3. But now, using the 

explicit results for the infinite strip of finite widthP^, we also show, as dashed curves in 

Fig. 3, the totally decoupled r = (or J2 = 0) plots for the two cases mi = 8 and 16. 

Clearly the uncoupled upper maxima fall below Tic just as do the coupled (r = 0.3) 

results. (It is worth remarking, however, that for a finite n x n Ising lattice with 

periodic boundary conditions, as studied by Ferdinand and Fishei^, the maxima in 

the specific heats lie above the bulk critical temperature Tic.) Nevertheless, there 

is clear evidence of a coupling or proximity effect in that the specific heats for the 

alternating, coupled system lie markedly above those for the decoupled (r = 0) strips. 

This same effect is seen in the experiments when finite boxes are coupled by a helium 
filnPZl. 

Complementary phenomena are observed around the lower maxima. Thus the 
dotted plots in Fig. 3 show the finite- width result for the situation r — )■ 00, or, more 
intuitively, Ji = 0, for nii = 8 and 16 (i.e., m2 = 16 and 32). These decoupled 
specific heats appear as very sharp, but still finitely rounded, spikes. However, it 
must be noted that these r — )■ 00 maxima lie below T2C, in accord with expectation 
for a finite-width strip. On the other hand, the maxima of the coupled alternating 
system lie above the limiting value T2C as seen clearly in the inset in Fig. 3. Once 



again there is an unmistakable proximity or enhancement effect that is found also in 
the experimental studies^. 

As a next step of our qualitative exploration, we present in Fig. 4 the effects 
of varying the relative separation s for significantly wide, mi = 18, strips spaced 
apart by weaker strips of relative strength r = 0.3 (as before). In this case the first 
point to notice is that Tc{r, s) increases quite rapidly towards Tic as the separation 
s approaches zero. Next, the uncoupled (r = 0) specific heats near Tic (shown 
dashed as in Fig. 3) all have maxima located at the same temperature, determined 
only by rrii = 18 for a finite width strips, while their magnitude is determined 
by s simply via normalization, either through relative area or on a per-site basis. 
However, there is still clear enhancement in the coupled layers even though the 
corresponding rounded maxima deviate very little in location from the uncoupled 
case. By contrast near T2c, as illustrated by the inset, the displacements even of the 
uncoupled maxima (shown dotted), depend significantly on the relative separation 
ratio s. Again, nonetheless, there is proximity induced enhancement of the peaks 
both in magnitude and displacement above T2c. 

Finally, we may enquire about the level of the specific heats around the bulk crit- 
ical point Tc or in the vicinity of the minima observed in Figs. 2-4 that lie roughly 
at Tmin ^ |(^ic + T2c)- One may ask, for example, how well the levels are ap- 
proximated by appropriately weighted sums of the uncoupled peaks around Tic plus 
some, perhaps reversed contribution from T2c. For these purposes, however, we need 
to proceed more quantitatively. 

III. SCALING EXPLORATIONS 

We would like to relate the observations embodied in Figs. 2-4 to more general 
scaling concepts. To this end, recalP^^^^^^ that a bulk system with a critical tem- 
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FIG. 3. The specific heats of alternating layered Ising models with relative strength r =0.3 
(as in Fig. 2(a)) and relative separation s = 2 for mi = 8,16. The dashed plots denote 
the corresponding decoupled specific heats when J2 = (r = 0), while the dotted lines 
represent the specific heats when Ji = (or r — t- cxd). The inset displays the distinct 
maxima near T2c, dotted below but solid above. 
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FIG. 4. Specific heat plots for relative strength r = 0.3 (as in Figs. 2(a) and 3) but with 
relative separations s = 1/3, 1 and 5/3 when mi = 18. Again the decoupled r = behavior 
is seen in the dashed plots while for the opposite limit, r — t- oo, the plots are dotted. The 
inset shows various coupled and uncoupled maxima near T2c- The short vertical lines locate 
the bulk critical points, Tc{r,s), which decrease as s increases. 
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perature Tc may be characterized by a correlation length ^{T) which diverges on 
approach to criticality as 

e(T)^eo/|ir with t = (T/T,)-1^0, (7) 

where i/ is a characteristic critical exponent while .^o is a length of order the lattice 
spacing a, or molecular size, etc. For 2D Ising systems one ha j^ * ^ * ^^ * ^ u = 1, whereas 
for superfluid helium in three bulk dimensions z/ ~ 0.67^. Then in a system limited 
in size by a finite length L = ia, the scaling hypothesis asserts, in general terms, 
that when i and ^{T) are large enough, the rounding of critical point singularities is 
primarily controlled by the ratio y = L/C,{T). 

Consequently, for the finite-size behavior of the specific heat per site, which di- 
verges in bulk as |t|~" where a is typically small (or even negative), the basic scaling 
hypothesis may be expressed as 

Cii;T)^r/^[Qix)-Qo]/a, (8) 

where Q{x) is the scaling function while the scaled temperature is 

X = t'^t oc y^'" = [L/e(T)]i/^ (9) 



and Qo > is a constant parameter. The exponent a in the denominator in (pj) 
allows for the limit a — )■ 0, which yields, with Q(0) — )■ Qo, a logarithmic singularity 
as is appropriate for 2D Ising systems. One may then take 

C{t, T) ^ [Qo/u) In £ + Q{^^''t) (10) 

as the basic hypothesis where, for use below, we note that at criticality one has 
C(£; Tc) ~ (Qo/t^) lni + Q{0). In fact, this hypothesis has been established explicitly 
for infinite Ising strips of width i and Q{x) has been explicitly determinecP^. 
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To apply these concepts to our layered Ising system, in the first case for the 
upper maxima near Ti^, we recall from Fig. 3 and Fig. 4 that leaving aside relatively 
small enhancements in magnitude, the total specific heat, C{Ji, J2',T), approaches 
rather well the limiting forms^ of a suitably normalized single strip of width mi. 
Accordingly, we subtract a contribution from non- coupled weaker Ising strips by 
defining 

Ci(Ji,J2;T) = (l + s)[C(Ji,J2;T)-C(0,J2;T)], (ll) 

where the normalization factor (1 + s) is needed for the scaling plots now to be 



examined. Finally in accord with ( 10 ) and the subsequent remark we introduce the 
upper or stronger net finite-size contribution 

ACi(Ji,J2;T) = Ci(Ji,J2;T)-Ci(Ji,J2;TiJ, (12) 

in which the value at the limiting critical point, Ti^, has been subtracted. If we 
accept the identifications i =^ mi and t =^ (T/Tic) — 1 and recall u = 1, we might 
expect ACi(T) to obey scaling in terms of the variable 

Xi = mi[{T/Ti,)-l]. (13) 

This expectation is very well supported by the plots in Fig. 5 for m,i = 18 and 
s = n/3 for n = 1, 2, ■ ■ ■ , 5: the "data collapse" is strikingly well realized. 

Beyond this, however, explicit calculations^ show that, asymptotically, ACi(T) 
is simply related to the limiting scaling function, Q°°{xi), for an infinite strip of 
coupling Ji and width rrii already known explicitly^. Specifically, allowing for nor- 
malization, yields 

ACi{Ji,J2;T)^{l + s)[C{Ji,0;T)-C{Ji,0;Ti,)]^Q^{xi)-Q'^{0). (14) 

Note that in this limit not only has the dependence on m,2 dropped out but also the 
dependance on J2. However, as regards the enhancement seen in Figs 2 to 4, we 
know that m,2 and J2 do play a role. This will be studied further below. 
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Let us now shift attention to the behavior of the specific heat peaks of the alter- 
nating system, near the lower (or weaker) limiting critical point, T2c- The rounded 
maxima are shown in detail in the insets of Figs. 3 and 4. Now we can follow the 



procedure that led to the definition (11). Thus we consider the normalized difference 

C2{Ji,J2;T) = {l + s-')[C{Ji,J2;T)-C{Ji,0;T)]. (15) 

Then, following again the previous analysis, the weaker net finite-size contribution 



may be defined as in (12), by 

AC2( Ji, J2; T) = C2( Ji, J2; T) - C2( Ji, J2; T2J. (16) 

It is natural to suppose that AC2(T) might obey scaling in terms of the new variable 

X2 = m2[{T/T2,)-l]. (17) 

This hypothesis is tested in Fig. 6 and, evidently, is remarkably successful, exhibiting 
excellent data collapse. But more remarkable yet is the evidence provided by the 
solid line plotted in Fig. 6. This derives directly from the limiting scaling function for 
an infinite strip of coupling J and finite width rrP^ but with the sign of the argument 



reversed. In other words, the previous asymptotic form (14) is now, as established 
irP replaced by 

AC2( Jl, J2; T) ^ (1 + S-')[C2{0, J2; f) - C2(0, J2; T2e)] ^ g°°(-X2) " g°°(0), (18) 

where the modified temperature T is simply attained by reflecting about T2c] explic- 
itly we have 

f{T) = T2c - (T - T2e) = 2T2, - T. (19) 

We may note, further, that in this limit the original dependence on both Ji 
and nil has vanished; but, once more, there are clear residual effects associated 
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FIG. 5. Plot of the upper or stronger net finite-size contribution, namely, ACi(Ji, J2; T) 
as defined in the text for strong strips of width mi = 18 at various separations but fixed 
r = 0.3. The solid curve is the plot of the specific heat of an infinite strip of width mi = 18 
and coupling Ji when its value at the bulk critical temperature Tic is subtractecP^. 
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FIG. 6. Plots versus the scaling variable X2 = m2[(T/T2c) — 1] of the lower or weaker 



net finite-size contribution AC2(Ji, J2; T), as defined in relations (|15[)-(16), for strips of 



width m,2 = 16 and r = 0.3 for various relative separations. The solid curve represents the 
corresponding asymptotic form Q°°(— X2) — Q°°(0), for an infinite strip of finite width and 
coupling J2 , with the temperature reflected about T = T2c~- The dashed curve represents 
data for a wider strip with ?7i2 = 60. 
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with the proximity and interlayer couphngs. To address these we may define an 
"enhancement" by subtracting from the total specific heat per site contributions 
deriving from the corresponding independent uncoupled strips. However, in doing 



this we must recognize — following Fig. 6 and the result (18) — that a reversed 



or refiected temperature variable is needed around T2c- To this end we utilize the 



modified temperature variable, T(T), via equation (19). Thus we define the net 
'' enhancemenf for fixed nii and m2 by 

S{Ju J2; mi, 7712-, T) = C{Ji, J2; T) - C( Ji, 0; T) - C(0, J2; f{r)). (20) 



It worth remarking parenthetically that in adopting this definition of the en- 
hancement we are, in particular, utilizing the theoretical results (I8p9) proved for 



the alternating Ising strips^. In more general situations (such as confined superfluid 



helium) the last term in ( 20 ) should be replaced by an symptotic term obtained by 
an initial data analysis of the behavior close to T2C such as led to the original (finite 
1712) form in Fig. 6. 

In Fig. 7, we plot the enhancement for our alternating Ising strips as a function of 
t = [{T/Tc) — 1] for nil = 8 and various separations s. One sees that the logarithmic 
divergence at t = is barely visible for s = 1, and essentially disappears for s > 1. 
In addition, as expected, the magnitude of the enhancement decreases as s (or 1712) 
increases; but by what law? 

To address this question we recall, first, that the leading correction to the asymp- 
totic form of the specific heat of an infinite strip of finite width m must arise from 
the two non- vanishing boundary free energy contribution^^^^'^^'^ which, yield a total 
specific heat term of relative order 1/m. The effect of this are already evident in 
Fig. 6 where the primary contribution (solid curve) is, especially for X2 > 2, more 
closely approached by the data for m2 = 60 than that for m2 = 16. It is clear that 
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FIG. 7. Plots of the enhancement £{t) versus t = (T/Tc) — 1 for mi = 8, r = 0.3 and 
various relative separations s. The short vertical lines above Tc{r, s), i.e., for t > are the 
corresponding positions of Tic, while below Tc{r, s), they locate T2c- 



such corrections must arise also in the bulk alternating strip system from the regu- 
larly spaced modified boundary or seams. By the same token, boundaries of surface 
effects should play same role in the experiments on the specific heats of small helium 
]^Q-jjg^il4i5i7i8i coupled via helium films. 

Accordingly, Fig. 8 presents the enhancements S{t) versus t oc [T — Tc{r,s)], but 



now multiplied by the factors (rrii + 77^2), which clearly corresponds simply to the 
density of seams in the bulk. It is striking that the maxima (close to Ti^) and 
the minima (near T2c) appear to rapidly approach constant values. This represents 
strong evidence that the enhancenment S{Ji, J2', mi, 7712] T) is of order l/(?7ii + 1712) 
as the relative separation, s = 7712 /mi, increases at fixed mi. 

However, by comparing Figs. 8(a) and 8(b), it becomes clear that the behavior 
of the rescaled enhancement peaks that approach Tic, when s increases, depend 
quite noticeably on mi, the width of the strong strips. Specifically, the enhancement 
peaks become both narrower, as indeed implied by Fig. 6, and taller as mi grows. 
Consequently, we will separately investigate the behavior of the enhancement close to 
Tic, noting that some logarithmic dependence on mi might be present. Nevertheless, 
for T < Tic Fig- 8 suggests that the enhancement rescaled by (mi + 7712) might 
approach a more of less constant shape when plotted against a suitable reduced 
intermediate temperature variable, say t^(T; Ji, J2; mi,m2). 

Then, since the expected scaling behavior must switch in the region between Tic 
and T2C we consider a definition that may differ somewhat according to the sign 
of t = [T/Tc{r,s) — 1]. To that end, we may recall the known values of the bulk 
correlation lengths, ^i(T) and ^2{T), for the strong and weak strips at the bulk 
layered critical point Tc{r, s), and note the {T = Tc) asymptotic identity (estabhshed 
iiP) that states 

Ci(Ji, J2;mi,m2) := (1 + s"^)[(Tic/Tc) - 1] 

^ ii\Tc) + ^2\Tc) := (1 + s)[l - {T2c/Tc)] = C2. (21) 

This suggests distinct but asymptotically equal rescaling factors for t above and 
below Tc{r,s). Finally, it is convenient to arrange that t"'"(T2c) = and t^{Tic) = 1: 
this yields the definition 

tt(r) =t/ci + (l + s)-i for t>0. 



19 



(m^ + ^2)^ 




FIG. 8. The enhancement £{t): (a) for mi = 8 as in Fig. 7, and (b) for rrii = 16, but 
muhiphed by {mi + m,2)- The short vertical hnes locate the corresponding upper limiting 
critical points, Tic- 
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t/c2 + (1 + sY^ for t < 0, (22) 

for which, as readily checked ,the interval t'''(Tic) — t^(T2c) is unity for all values of 
Ji, J2, fni and m2. 

Using this intermediate reduced temperature variable, Fig. 9 presents the rescaled 
enhancement in the range below t"'' c^ 1.5 down to t"!" ~ —0.5, where S{T) is negative 
and below which little remains to be seen. Note, however, that the rounded peaks 
close to t^ = 1 have been omitted. While some peculiar rapidly varying features 
appear in the interval for t"'' = to 0.1 associated with the limiting behavior of the 



subtracted third term in the enhancement definition (20), the overall set of data 
support the hypothesis of a unversal limiting plot when s — > 00 for large mi. 

Finally, we investigate the enhancement in the vicinity of Tic, recalling first, the 
subtracted and scaled total specific heat plots presented in Fig. 5. The maxima of 
the enhancement as mi runs from 8 to 64 are studied in Table I. Table I lists the 
values of (mi + m.2)Sic{m.i) evaluated at Tic and (mi + m.2)£max{jfii) together with 
their ratio. Extrapolation of the ratio suggests a finite limit £max/£ic — 1-06. This 
is consistent with a finite-size scaling expression for the asymptotic shape of the 
enhancement peaks. However, it is clear that the values of {mi + m2)£ic rise steadily 



with mi. As already hinted, in light of the leading asymptotic behavior relation (10) 
a logarithmic factor in the correction terms must be anticipated (although seemingly 
absent for T < T2C as evidenced by Figs. 8 and 9). In fact the slightly extrapolated 
fit 

(mi + m2)^ic(mi) ~ 0.0800 ln(mi/mo), m^ = 0.36, (23) 

seems optimal as may be judged from the last column of Table I, which records the 
ratio TZ{mi) of the two sides of this expression. 

Last, to study the asymptotic shape of the enhancement peaks, note from Fig. 
5 that the associated total specific heat peaks near Ti^ are displaced by a known 
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FIG. 9. Plots of the rescaled enhancement (mi + ra2)£{t) for mi = 8, 16 and 32 and 
separations = 1/2,1,2,4,8 versus the reduced temperature variable, V{t), defined in the 
text. The color code is black for mi = 8, blue for mi = 16 and red for m,i = 32. 



scaled amount. Specifically, the scaled asymptotic maximum located at z™*^^ = x\ = 
—0.7869^. The corresponding values of xi for the maxima of the enhancement as 
nil increases are listed in Table I together with the ratio x^""^ {nii) / x\. The data 
seem subject to numerical imprecision but inspection of the behavior of the ratio, 
say vs. l/rrii, we may reasonably conclude that the enhancement peaks, while of 



22 



somewhat different shape, with x'^"'^{mi)/xi — )■ 0.6, nonetheless, obey asymptotic 
finite-size scahng in terms of the strong bulk correlation length C,i{T). 



mi 


(mi + m2)£ic 


(mi + m2)£max 


^max / ^Ic 


n{mi) 


^max 


^max/3.0 


8 


0.24884 


0.34034 


1.368 


0.9989 


-0.805 


1.02 


16 


0.30386 


0.36827 


1.212 


1.001 


-0.917 


1.17 


24 


0.33628 


0.39241 


1.167 


1.001 


-0.755 


0.959 


32 


0.35936 


0.41036 


1.142 


1.000 


-0.662 


0.842 


40 


0.37731 


0.42505 


1.127 


1.000 


-0.741 


0.941 


48 


0.39197 


0.43756 


1.116 


0.999(9) 


-0.731 


0.929 


56 


0.4044 


0.4486 


1.109 


0.999(7) 


-0.698 


0.887 


64 


0.4152 


0.4579 


1.103 


0.999(5) 


-0.62 


0.79 
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